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A microscopic Hamiltonian theory of the FQHE developed by Shankar and the present author 
based on the fermionic Chern-Simons approach has recently been quite successful in calculating 
gaps and finite tempertature properties in Fractional Quantum Hall states. Initially proposed as a 
small-g theory, it was subsequently extended by Shankar to form an algebraically consistent theory 
for all q in the lowest Landau level. Such a theory is amenable to a conserving approximation 
in which the constraints have vanishing correlators and decouple from physical response functions. 
Properties of the incompressible fractions are explored in this conserving approximation, including 
the magnetoexciton dispersions and the evolution of the small-g structure factor as i/ — > i . Finally, 
a formalism capable of dealing with a nonuniform ground state charge density is developed and used 
to show how the correct fractional value of the quasiparticle charge emerges from the theory. 
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I. INTRODUCTION 

The Fractional Quantum Hall (FQH) effectB has in- 
trodutied us to new states of electrons in high magnetic 
fieldsH. Laughlind explained the essence of the effect for 
fractions v = l/(2p+ 1) and showed that the electrons 
form strongly correlated states which are incompressiblie 
and have quasiparticle excitations with fractional chargecl 
e* =e/(2p+l). 

A unified understanding of all fractions v = p/ {2sp 
was achieved by the Composite Fermion picture of Jair 
In this picture, the true quasiparticles are electrons 
dressed by 2s units of statistical flux, which are called 
Composite Fermions (CFs). At a mean field level, the 
CFs then see a reduced field B* — B/ {2sp+ 1), in which 
they fill p CF-Landau levels (CF-LLs), and exhibit the 
integer QHE. This picture has been very successful in 
describing the FQH statea3. 

Numerical diagonalizatiQiis of small systemajfl or the 
use of trial wavefunctionsB'S are the best approaches for 
calculating many properties of FQH systems such as 
ground state energies, gaps, magnetoexciton dispersions, 
ground state spin transitions, etc. However, there are 
other interesting properties which are not accessible to 
these approaches, the most important of which are dy- 
namical response functions. Partly to address these is- 
sues, and also to gain a better understanding of FQH 
states, several field-theoretic approaches were developed 
during the past fifteen years. Most of these approaches 
are based on the Chern-Simons(CS) transformationci, 
which performs flux atjtaciiHient via thCr-CS gauge field 
to obtain either bosonaij^E^ or fermionstil. These theo- 
ries have provided us with a link between the microscopic 
formulation of the problem and expepLipfint, both for in- 
compressible and compressible statesEjO. 

Mpt|ivated by some outstanding issues in previous 
workO, Shankar and the present author recently deveL- 
oped a hamiltonian CS theory for the FQH statesE^Ej. 



Inspired by the work of Bohm and Pine on the 3D 
electron gas, we enlarged the Hilbert space to introduce 
n high-energy magnetoplasmons degrees of freedom, (n 
also being the number of electrons) at the same time 
imposing an equal number of constraints on physical 
states. Upon ignoring the coupling between the oscil- 
lators and thc,jfermions we obtained some well known 
wavefunctionsBu. A final canonical transformation was 
employed to decouple the fermions from the oscillators 
in the infrared limit. 

After the final canonical transformation, we obtained 
the following expressionstj for the electronic charge den- 
sity pe and the constraint x for = p/2p + I at small 

q- 
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where A, refer to the annihilation and creation oper- 
ators of the magnetoplasmon oscillators, I = 1/VeB is 
the magnetic length, Hj = Pj + eA*{rj) is the velocity 
operator of the CFs, and c — \phj. Instead of using the 
expression of eq(|^) for the electronic charge density we 
added a multiple of the constraint to it in subsequent 
calculations. The rationale for doing so is that in an 
exact calculation in the physical subspace (which is an- 
nihilated by the constraints) using the operator p^. — 
with arbitrary ^ instead of pe should not change the phys- 
ical answer. However, in an approximate calculation in 
the full space different choices of ^ will produce different 
answers. The choice ^ = = 2v is particularly advan- 
tageous, as will be seen below. With this "preferred" 
choice, the expression for the density becomes 
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The oscillator piece saturates Kohn's theore m0. The 
rest, to be called pp, has the following properties: 

• The presence of the velocity operator 11 in pp shows 
that the CFs respond to the effective field B* . 

• Pp satisifies the magnetic translation algebra 
(MTA)E3 to lowest leading order. Since this is the 
algebra of the electron density in the lowest Landau 
level (LLL), we have performed the LLL projection 
correctly in the infrared. 

• Note that pp is a sum of a monopole with charge 
e* = e/ (2p+ 1), which is the-charge associated with 
the CF, and a dipole pieceES which alone survives 
aiv — 1/2. (A number of recent-coajStructions have 
emphasized this dipolar aspectcJ^EZl) . 

• We also find that that as g ^ all transition matrix 
elements of pp from the HF ground state vanish at 
least as (a conseqi^ce of Kohn's theoremO for 
incompressible statesEj) . 

These are all nonperturbatively correct properties of 
the true electronic density operator in the LLL. In our 
formulation the full Hilbert space of the CFs is too large, 
and the physical LLL will be reached only after the con- 
straints have been properly imposed. Unfortunately, we 
do not know how to select the physical states opera- 
tionally in the final representation. The use of pp rep- 
resents a shortcut that obviates the need to impose con- 
straints. The idea is that since pp displays the same prop- 
erties in the full space that the true electronic density is 
known to possess in the physical space, calculations car- 
ried out with pp in the full space will hopefully provide a 
reasonable approximation to the true physical answers. 

The Hamiltonian of the low-energy sector now is 

H=\j ■0^v{q)-pp{-q)-pp{q). (4) 

where v{q) is the electron-electron interaction. 

Once the constraints have been used to construct pp, 
they are ignored in subsequent calculations. One assumes 
that the Hartree-Fock(HF) approximation is sufficient to 
capture the correct physics-j-This assumption is borne 
out by calculations of gapsEaE^ aad-the temperature de- 
pendence of the spin polarizationE3tII. It even provides a 
reasonable iqiialitative description of the magnetoexciton 
dispersionsS 

Despite these successes, the above formulation has a 
few disadvantages. Firstly, the expressions for the den- 
sity and constraint operators are known only to first or- 
der in q. The use of the small-g preferred density oper- 
ator for all q renders the gap infinite for the Coulomb 
interaction, and one has to use a potential that is ex- 
ponentiaUjj suppressed at large q in order to get finite 
resulta^J'Eil. In computing magnetoexciton dispersions 



using the small-g approachE3 the author used a cutoff in 
the number of CF-Landau levels (CF-LLs) corresponding 
to keeping the correct number of single-particle states in 
the Hilbert space to obtain finite results for the magen- 
toexciton dispersions for any interaction. However, for 
realistic potentials the disper|Sipns were still strongly de- 
pendent on the CF-LL cutoffE3, casting doubt on the ro- 
bustness of the approach. A related disadvantage is that 
the density operator satisfies the MTA only to leading 
order in q. 

As stated above, the preferred density oper|ator was 
obtained by adding a multiple of the constraintllj to the 
canonical transform of the original electron density. This 
made certain nonperturbative features of the charge den- 
sity explicit. However, it would be desirable and more 
satisfying to see these nonperturbative features emerging 
from the famralism rather than having to put them in. 

Shankart3 extended the previous small-g approach to 
all orders in q by exponentiating the first two terms of 
the small-g expressions for the density (eq(|l])) and the 
constraint (eq(||)). The resulting expressions (considered 
in more detail in the next section) have some very de- 
sirable algebraic properties; the density operator obeys 
the magnetic translation algebra exactly at all q, and the 
density and the constraint commute. This implies that 
the Hamiltonian constructed from the density operator 
commutes with the constpdnt, which makes a standard 
conserving approximationE3 possible. In such an approx- 
imation, the correlator of the constraint with either the 
physical density or itself is explicitly zero. It is worth 
noting that the algebra found biz_Shankar reduces to that 
found by Pasquier and HaldaneEj for bosons aX, v = 1 (in 
this case = v m. the above formulas, and bosons at 
V — 1 are mapped to fermions in zero average field, just 
as electrons are at v_= i). This case was treated in 
great detail by ReadEj, who computed bosonic response 
functions in a conserving approximation. 

In this paper, we will follow this conserving route 
to calculating response functions for the incompressible 
fractions. This corresponds to using the HF approxi- 
mation on the Hamiltonian to calculate single-particle 
energies and calculating bosonic correlators in the time- 
dependent HF (TDHF) approximation. Several results of 
interest are obtained. Firstly, it is known that the spec- 
tral function of the correlator of the projected electron 
density operator cannot contain terms of order q^, since 
this would violate Kohn's TheorcmEJ. Naively, a den- 
sity perturbation produces excitons from filled states to 
empty states. The electronic density operator has ma- 
trix elements of order q to some of these states, and 
it is interesting to see how these disappear in the final 
correlator. It turns out that there is always one linear 
combination of naive excitons with a zero eigenvalue in 
the TDHF approximation. This corresponds in a very 
close sense (to be made precise later) to the constraint, 
and represents the unphysical sector of the theory. The 
electronic density operator does not couple to this un- 
physical sector in the TDHF approximation, and does 
indeed have matrix elements of order q^ or higher to the 
physical excitonic states, in compliance with Kohn's the- 
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orem. As a consequence of the matrix elements the 
structure factor S{q)^ which is the equal time correlator 
< pe{q,t)pe{—q,t) >, should go asu^'' in the incompress- 
ible states. However, Read's workE 



I shows that at ly — ^ 
this behaves as q'^logq. We will show that the for the 
fraction v = p/ {2p +1), the coefficient of the g'* term in 
the structure factor diverges (for p > 1) as 
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is one of the 



This result shows that as p — > oo, a power series expan- 
sion of S{q) must-fail, and is consistent with the above 
nonanalytic resultnj. 

We will also present results for the magnetoexciton 
dispersions for 1/3 and 2/5. The positions of the min- 
ima are found, to be exactly at the values of ql^ found 
numcricallyBia. This attests to the fact that although 
the all-g formalism was derived by extending a small-g 
approach, it nevertheless manages to capture the struc- 
ture of the magnetoexcitons up to ql^usi- 

The fractional quasiparticle chargajO' 
most striking and robust features of the FQHE, and it 
must emerge from any theory that correctly describes 
the FQHE. To see how this appears in the all-q Hamilto- 
nian theory, a further extention oLthe formalism is nec- 
essary. Shankar's all-g formalismE3 is designed to work 
best in a uniform background. In an inhomogeneous 
background there will be a nonuniform average of the 
electronic charge density n{r) = n + Sn{r), which we ex- 
pect by the principles of flux attachment to lead to a 
nonuniform effective magnetic field. Clearly, adding a 
localized quasiparticle produces just such a nonuniform 
charge density and effective field. A version of the all-q 
theory applicable to first order in Sn/fi will be presented, 
and used to show the emergence of the fractional quasi- 
particle charge. Such a formalism is also useful if one 
wants to treat disorder or edges in a conserving approx- 
imation. 

The structure of this paper is as follows: In Section || 
we provide a brief introduction to the Hamiltonian the- 
ory of the FQHE, including the single-particle states. In 
Section III we will present the structure of the conserving 
approximation and the details of how it is carried out. 
In Section IV we will analyse the TDHF matrix, show 



that it has zero modes corresponding to the unphysical 
modes. The result for the small-g structure factor for 
ly — p/{2p + 1) will be derived in Section 0. In Section 
VI we will present results for the magnetoexciton disper- 



sions for i 



and f for q not too small. In Section VH 



we 



will present the extention of the all-q formalism to the 
case of nonuniform density, and see how the fractional 
quasiparticle charge emerges naturally in the theory. We 
end with some conclusions. 



II. HAMILTONIAN THEORY OF THE FQHE 

As described in the int rodu ction, the small-g version 
of the Hamihonian theoryllalla of the FQHE was derived 



by introducing extra degrees of freedom which assumed 
the role of the magnetoplasmons, and supplementing the 
Hamitonian with constraints. After decoupling this led 
to a theory with many desirable properties, but also some 
deficiencies, notably involving the consistency of the the- 
ory for all q. 

We will take the alLw extention of this formalism as 
described by Shankaitj as our starting point. In this 
formalism two sets of guiding center coordinates are as- 
sociated with each CF. 
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Here = 2t^, and r and H are the position and velocity 
operators of the CF. The algebra of these objects is of 
fundamental importance. 



[Re,a, Re.fs] — —il'^^aji 
[Re,a,Rv,p\ = 

[Rv,a,Rv,p\= i-^£at3 



(8) 
(9) 
(10) 



The first set of coordinates Re is identified with the 
guiding center coordinates of the electron, and the sec- 
ond set Ry corresponds to the guiding center coordinates 
of an object with a magnetic algebra charge of — e/c^. We 
will call this object the psaudmjortex (for brevity, Shankar 
simply calls it the vortexEjEJ). Upon projection to the 
physical subspace we expect that this coordinate will dis- 
appear, leaving behind the proper correlations between 
electrons, that is, vortices. Magnetic translation opera- 
tors Tej Ty and densities Pe, py can be formed from the 
two sets 

i 

py(q) - e"i"'"/^''"Ty{q) = e-?'''/^-' ^ g-'^-^-' (12) 

i 

The operator Tf. obeys the algebra of the electronic 
magnetic translation operators projected to the LLLcj 

[T,{qi), T,{q2)] = 2i sin(^gi x q2)TM + 9*2) (13) 
while the operator Ty satisfies a similar algebra 

[Ty{qi),Ty{q2)] = -2isin(^gi x q2)Ty{qi + q^) (14) 

The density formed from the electronic guiding center 
is identified with the electronic density operator, while 
that formed, from the pseudovortex guiding center is the 
constraintEa. It is easy to verify that the small-g limit 
of these expressions indeed reduce to those derived in 
the small-(7 theory (eqs(|l|,||)). Since they commute with 
each other, the Hamiltonian formed from pe alone will 
commute with the constraint. _ 

In previous uses of this all-q formalismEil, the preferred 
density pp = pe — c^Py was used. The preferred den- 
sity has many desirable features, as described in the in- 
troduction, and gives numbers in good agreement with 
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numerical, and, fixperimental results for many physical 
propertiesE3E30. However, in addressing issues of prin- 
ciple, such p^Ae form of the density correlator at small-q 
for V ~ 1/2E3c3, naive calculations based on the preferred 
density fail to produce the correct answer, and a conserv- 
ing calculationESES is needed. 

Before we can compute density correlators, we need 
to construct the single-particle states of the theory. Our 
Hamiltonian is 



;v{q)pe{q)Pehq) 



(15) 



From eq(|6|) we see that the coordinate Re contains H, 
the velocity operator in the effective field. Thus, Landau 
level states in the effective field make a natural basis in 
which to represent the various operators. In the Landau 
gauge, these states are 

= Ce^^y/(n\-i^-xfmrrHU^^ _ x)/i*) (le) 

where C is a normalization constant, are the Her- 
mite polynomials, /* — \/\/eB* is the magnetic length in 
the effective field, n represents the Composite Fermion- 
Landau Level (CF-LL) index, and X — 2nj{rY/L is 
the degeneracy index (the degeneracy of each CF-LL is 
2^(;.)a ). In terms of these states the translation opera- 
tors can be represented as follows: 



Te{q) 



E e 



Qyl* 



n2,X + 



Qy' PKira2 (q) (17) 

2 

Qy^* Xnin2 io) (18) 



where the matrix elements are 



-n2)(e,-;r/2) 



pn.«.(9i = (-ir<+"% 

(cQ/V2)l"i-"^le-^''3'/4iL"r"''(c'QV2) 

X^{q) = (-l)"<+"2y^g»(ni-n2)(0<,-7r/2) 

(Q/cV2)l"i-"^le-Q'/4-'LL"<^-"^l(QV2c2) 



(19) 



(20) 



Here c — \/2i^, Q = ql* , n< (n^) is the lesser (greater) 
of ni, n2, and Ljj are the Laguerre polynomials. 

It can be shown that the CF-LL states qr'?. ^ct 
Hartree-Fock (HF) states of the HamiltonianHEJ, that 
is, they arc the best possible Slater determinantal states 
for this problem. The single-particle HF energies are 

1 r (fq 



= o / 7^v{q)e-''''^'^Y.^l-'2NF{m))\pr,ra{q)? 



(21) 

Note that the energy depends only on the CF-LL index, 
and not on the degeneracy index. The first term cor- 
responds to the Hartree energy, while the second term 
is the Fock energy. When the preferred density is used 
the Hartree contribution is nontrivial, and depends on n. 
However, when the electron density of eq(|l|) is used, as 
we do here, the Hartree energy is a constant independent 
of n (the methods of the Appendix can be adapted to 
show this). 



III. THE CONSERVING APPROXIMATION 

To illustrate the idea of a conserving approximationll, 
let us go back to the Hamiltonian constructed from the 
all-q density operator Pe{q), eq([l5|) . This should be sup- 
plemented by constraints that demand that the pseu- 
dovortex density should annihilate all physical states 



Pviq)\^phys >= 



(22) 



Now define a time-ordered Green's function Gy^iq^t) 
as follows (one can similarly define Gee and Gm) 

G^e{q,t- t') = -i < Tp^{q,t)pe{-q,t') > 

= -iOit-t') < p,{q,t)pe{-q,f)> 

-iQ{t' -t) < pe{-q,t')p,{q,t)> (23) 

where T is the time ordering operator. Consider the time 
development of this Green's function 

-i^G^e{q,t-t') = -5{t-t') < [p,{q,t),pe{-q,t')] > 
-i<T[H,p,{q,t)]pe{-q,t') > 

Since p^ and pe commute, p^ also commutes with H. One 
immediately sees that G-ue is a constant. If one sets the 
constraint to zero initially, then it remains zero, and all 
its correlators also remain zero. 

The above is true as an exact statement. However, in 
practice one uses some approximation scheme to calcu- 
late the Green's function. In our case, an approximation 
that respects Gye = Gw = will be conserving. Let us 
see what a natural approximation scheme might be. Con- 
sider the time-ordered Green's function Gee {q, t) defined 



Gee{q,t-t') = ~i< Tpe {q,t)pe{-q,t') > 

and look at its time development. We have 



d 

l-Q^Gee{q,t) 



S{t) < [peiq,t),pei-q,0)] > 
-i<T[H,peiq,t)]pei^q,0)> 



a2"/3l ^ "ai"/32 



(24) 



(25) 



(26) 



From the commutation relations of pe with itself (eq([l3|)), 
it can be seen that a Green's function involving three 
densities will arise. Further time derivatives will involve 
higher-order density correlators. The natural way to 
truncate this hierarchy is to make a mean-field approx- 
imation at some stage that reduces a product of den- 
sities to a single, density. One of the simplest of such 
approximationscj reduces [H, pe] , which is a product of 
four fermi operators, to a product of only two, by using 
the averages 

- < di.dp^ > dl^df3, - < dLda, > dt, d 



(27) 



Here < dl^dp >— Saf:iNp{a), where Np{a) is the Fermi 
occupation of the single-particle state a. 
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Using the HF states and their occupations in the above 
truncation is known as the time-dependent HF (TDHF) 
approximation. In the electron gas it is known that this 
approximation is conserving (in that case the consecra- 
tion laws to be satisfied result from gauge invarianceE3) . 
ReadEj has shown that this approximation is conserving 
for the V — \ boson problem as well. We will explicitly 
show below that TDHF is conserving for all principal 
fractions. 

Let us proceed to some more explicit details. 



A. TDHF Formalism 

The final all-g theory has mapped the original strongly 
correlated electron problem into a weakly correlated CF 
problem with constraints. This CF problem has LL struc- 
ture in the effective field, and integer filling for CFs. 
All the techniques that have been brought to bear on 
the problem of computing magnetoexcitons in the inte- 
ger quantum Hall (IQH) states can be applied to this 
problem. The computation of magnetoexciton disper- 
sions has a long_|historyc3 in the IQHE. Many different 
approximationaljcd have been employed to obtain these 
dispersions in the IQHE, most jOf jSSihich can be subsumed 
into a unified TDHF treatmentll3H. In this approach one 
allows the hamiltonian to scatter particle-hole excitations 
with different Landau-level indices into each other, and 
diagonalizes the resulting matrix. As mentioned before, 
this method has been applied-, to the FQHE using the 
small-(? hamiltonian approached. Most of the details can 
be found there, so the following will touch on the high- 
lights of the formalism. 

The TDHF matrix can be derived from many differ- 
ent points of view. We will find it most convenient r-to 
follow the equation of motion approach of RickayzenE3. 
We start with the magnetoexciton operator 

= ^ (28) 

^ — ' ra\.X ^ — m2,v\-| ^ — 

X ^ 

This represents a CF-particle-hole excitation with a 
definite momentum q. To avoid a proliferation of indices 
we will use the shorthand — m\m2 (and sometimes 
V = nin2, not to be confused with the filling factor!). 
Now we define a time-ordered Green's function of two of 
these operators 

Gi^I■,^i';q■,t-t') = -i< TO^iq,t)0^,i-q,t') > (29) 

Taking the time derivative we obtain 



f G(m;a*';?;0 - < T[H,o^{q,t)]o^,{-q,o) > 

-5{t) < [O^{q,t),O^,{~q,0)] > 



The average value of the equal time commutator in the 
ground state (the final term of eq(^0|)) can be carried out 
to give 

{Npimi) - NfM) (32) 

Now the Hamiltonian is expressed as a four-fermi opera- 
tor 



q,{X,}{n,} 



I'd 



(33) 



One performs the commutator O^] , and in the TDHF 
approximation one contracts the resulting four-fermi op- 
erator to give a two-fermi operator. After the sums over 
degeneracy indices are performed one finds 

[H, = (e(mi) - e(m2))0^(-z) + 

{Np{m2) - Np{m,)) ^ 0.{q) (^-^e-^'^^' p^.r.. (q) x 

p™,™, i-q) - J ^v{s)e-"'"/^pr.,rn, (s)p™.„. {^s)e^^^'^ ] (34) 

One can now define a vector ^'(/i), which corresponds to 
the operator 



(35) 



Thus the action of commuting with H on in the 
TDHF approximation can be represented as the right- 
multiplication of ^' of a matrix 7i(/i; v] q) in the space of 
these vectors 

TL{n\v]q) = (5„ii„i5„i2„2(e(mi) - e(m2)) -I- 

I ^vis)e-^^'^'/'p^,„^A^)p^.nA-^e^''"^'^ (36) 

It is clear that finding the eigenvalues and eigenvectors 
of this matrix will also enable us to solve for the Green's 
function. Assume that one has found the right and left 
eigenvectors and corresponding eigenvalues, labelled by 
a 



H{p; V, q)-^'^{v; q) = E^{q)-^^{y- q) 
'f^ip; q)n{p, V- q) - E^{q)^^Ju- q) 



(37) 
(38) 



(30) 



Now one Fourier transforms G to convert it to a function 
of frequency lo according to 



where sums over repeated indices are implicit. 

We can choose to normalize the eigenvectors in the 
conventional way 



oo 

G{t)= J ^e-'G(c.) 



^^ip;q)^^{p;q)^d^f3 



(39) 



(31) 



If the matrix has a complete set of eigenvectors (we 
will see in the next subsection that this is not true, and 
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see how to deal with this situation) then one can expand 
(w - n)-^ as 

to — hin 



Now the solution to the Green's function can be ex- 
pressed as 

G(^l■,^i'\q]UJ) = 2;^E 

*«(/^;9l;:i*^(™Xi;'?)(^F(™i) "iVi.(™i)) (41) 



IV. FORMAL ANALYSIS OF THE TDHF 
MATRIX 

In this section we will first show that there is a definite 
connection between the right and left eigenvectors, and 
that the eigenvalues of H. always come in pairs ±Ea- 
Next we will show that there is always one eigenvector of 
TL with a zero eigenvalue, and that this corresponds to the 
constraint. There is also another zero eigenvalue which 
does not correspond to any eigenvector, and the matrix 
H. is defective, that is, does not have a complete set of 
states. We will see that we can nevertheless compute all 
the physical Green's functions, since they do not couple 
to the unphysical sector. 



A. Two Useful Properties of Ti. 

It is easy to see that TL is neither a symmetric not a 
Hermitian matrix, so in general one should expect right 
and left eigenvectors which are not Hermitian adjoints 
of each other. However, there is a definite connection 
between the right and left eigenvectors with the same 
eigenvalue: 

^^{v;q) (X [NFin^) ~ NF{n,)){^^M,q)y (42) 

This can easily be verified by noting that Pnin2{<l) — 
{Pn2ni{—<l))* 1 and taking the complex conjugate of the 
eigenvalue equation for the left eigenvector (eq(p8[)). In 
fact, any multiple of the above will form a good right 
eigenvector. However, to satisfy the normalization con- 
dition eqp9), we need to take 



(43) 



The second property is that eigenvalues occur in equal 
and opposite pairs ±Ea ■ By focusing on the direct term, 
it is easy to verify that the phase of a term in Ti is 

T-ld{fJ.;i^;q) _ / i \ri<+m<+mi+Tii 
g^(e,-7^/2)(nl-Hm2-n2-ml)(7VJ.(m2) " A^F(mi)) (44) 



It is clear from eq(|36|) that fi such that iVp'(m2) = 
Npinii) are not coupled to any other modes, and there- 
fore have a time-evolution given solely by their energy 
e(mi) — e{ni2). The nontrivial part of Ti consists of 
those vectors for which A''^? (777,2) — ^f("7i) 7^ 0. We 
can divide this subspace into a positive energy part with 
NF{m2) — Npimi) > 0, and a negative energy part with 
NF{m2) — NF{mi) < 0. Clearly, for the trivial part 
of the spectrum, one can immediately construct a nega- 
tive energy mode Om2mi for every positive energy mode 

17712 • 

Before we proceed, let us note that rotational invari- 
ance guarantees that the eigenvalues Ea{q) depend only 
on the magnitude of q. In order to analyse the structure 
of Ti the choice 0g = is particularly convenient. With 
this choice, it is easy to see that 



Tld{mim2\nin2]qx 



-TCd{m2mi; n2ni; q^ 



(45) 



With a little more effort it is possible to show that the 
same relation holds for the exchange part, and thus for 
the entire matrix Ti, which has the following structure 



TL^ 



TL++ TL+- 
-TL+- -TL++ 



It follows that if 



(u+ U-)TL = E{u+ u_) 



(46) 



(47) 



is a left eigenvector of TL with eigenvalue E, then 

(u_ u+)TL = -E{u- u+) (48) 

is a left eigenvector with eigenvalue —E. The only excep- 
tion that could occur is if = 0, in which case (u_ U-^) 
might possibly be the same as {u+ U-), and one is left 
with only a single eigenvector. This case is realized in 
our conserving approximation, and we now turn to its 
analysis. 

B. Unphysical States in TI 

The most important property of Ti is that it always has 
at least one left eigenvector with zero eigenvalue, namely 



(77717772; g) = X7nim2(9) 



(49) 



The requisite formalism to show this explicitly is de- 
tai led i n Appendix I, from which we borrow two identities 
eq(p6l) 



^Xmim{qi)Pmm2{q2) = 



< 777i|e"''(^+'=«^)--^CF|TO2 > e-^'^ixQa 



(50) 



and eq(137) 



^Xmim{qi)Pmm2{q2) =e ^ p,77i7n (^2 )X7n7772 (^l ) 



(51) 
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Let us now consider 



mim2 



27r(Z*)2 



{q)H{mim2;nin2;q) = (e(ni) 

77117712 



e(n2))X«in2('Z) 



77117772 



XtTII 7712 P711712 



(52) 



Let us consider the direct and exchange terms separately. 
In the direct term, one m index can always be summed 
freely, while the other is constrained by the Fermi occupa- 
tion factor Np. The sum over the free m gives, according 
to eq(|0|) 

^NF{m2)pn,nAq) < ^2 |e-'( ^ ^Ci^ |m2 > 

7112 

-J2NF{mi)pr,,nM < mi|e"'(i-=)?-^ci.|rni > (53) 

mi 

The two terms are immediately seen to cancel. Now let 
us turn to the exchange terms, and consider the one that 
has the factor jVp (777.2), a-nd a free sum over 7771. In this 
term, one can use eq(p^) to exchange the p and x matrix 
elements to obtain 



{2n) 



Xnimi 



(54) 



Notice that the phase factor e*"^^*^ has been cancelled by 
an opposite phase factor from eq(|5l]). Now the angular 
s integral forces 7771 = n2 for a rotationally invariant 
potential, and the result contains the Fock energy of the 
state 772 



' {n2)Xnin2i(f) 



(55) 



Similarly, the other exchange term proportional to 
Np{mi) ends up giving — e^(77i)xnin2 (?)■ Due to the 
peculiar nature of the Hamiltonian the Hartree energy 
is a constant independent of the CF-LL index, and the 
difference of the Fock energies is the same as the differ- 
ence of the full HF energies. Thus, the exchange contri- 
butions cancel the diagonal term (e(r7i) — e(r72))xnin2 {q}^ 
and Xnin2 i'f) is indeed an left eigenvector with zero eigen- 
value for 7i, as claimed. 

If we try to obtain another eigenvector with zero eigen- 
value using eqs( ^ , ^8|) , we are foiled by the fact that 
Xnin2(fe) = Xn2"i(fe), which givcs us back the same 
eigenvector. The matrix Ti. is defective, that is, it does 
not have a complete set of eigenvectors. The null sub- 
space of Ti. is two-dimensional, one of the vectors being 
the eigenvector (77717772; = Xmim2('?)- The other 
vector (caU it '^q,) is such that *q,7^ = ^q- In other 
words, in this subspace, H looks like 



(56) 



Fortunately, this singular behavior of the matrix is not 
a problem in computing physical correlators, because the 
density does not couple to the null subspace, as will be 
seen explicitly in the small q limit in the next section. To 
compute physical correlators, one need only project out 
this two-dimensional null subspace, and then the matrix 
inverse (uj — H)~^ looks exactly as in eq(|4l|), except that 
the sum is only over physical eigenvalues. 



V. SMALL Q ANALYSIS OF H 

The goal of this section is to compute the small-g be- 
havior of the structure factor S{q), which is the equal- 
time density-density correlator in the ground state. To 
this end, we will carry out an analysis of Ti for small q, 
and show that one needs to keep at most a 6 x 6 matrix 
to get the structure factor correct to order g^. One of 
the results of this section will be that the coefficient of 
this q"* term diverges as p — > 00, showing that the limits 
q ^ and p ^ 00 do not commute. 

Noting that the density can be expressed as 



Peiq) 



77117712 



PtTI 17112 7171 7172 



(57) 



we can write the density-density correlator S{q,uj) =< 
Pe{q,uj)pe{-q,~uj) > as 



S{q,uj) 

Pm\m2 



a, 717 17112 

(771 772 : q){NF{n2 



{q)^^{f,;q)—. 



(58) 



' uj — Ea-^iesgn{Ea) 

NF{ni))pn2m{-q) 

where the sum is over physical states only, and the ie 
factors appropriate to the time-ordered Green's function 
have been inserted. Now one can go back to the time- 
domain Green's function, and take the limit of equal time 
to write in a transparent notation 



S{q) = 



E <p{q)\^Siq)><^t{q)\p{q)> (59) 
where N is the number of particles, and we have defined 

(^ P('Z)l)mim2 Prnirn2 

(q) (60) 
{\p{q) >)mim2 = (NF{m2) ~~ N F{mi)) p^,^^^ {-q) (61) 

Recalling the relation between the right and left eigen- 
vectors eq(^), and noting that N = pL'^/2Tr{l*)^ for 



— _E. 

2p+ 



j, we finally obtain 



Siq) = -e 
P 



J2 |<p(9)|*f(<7)> 



(62) 



Our task in this section will be to compute all terms in 
this sum that contribute to order {ql)'^. It turns out that 
the case of a single CF-LL filled is nongeneric, while any 
p > 1 — p/2p + 1) is generic. We address the simpler 
p = I case first. 
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A. The Structure factor for the Laughlin fractions 

From the structure of eq(|6^) it is clear that only mag- 
netoexciton operators with Np{mi) ^ Np{m2) need to 
be kept in computing the density correlator. For v = ^ 
this corresponds to Omo (positive energy) and Oom (neg- 
ative energy). We will assume that the eigenvectors 
are analytic at small g, which turns out to be correct in 
the incompressible states. Furthermore, we will need the 
matrix elements p,„o which are 



n{mid : nO : q) = 



v{q) 



Pmoiq) 



(63) 



-(-l)"<^™'"^+'"(sgn(n - m))"-" / ^v{s)e-s"/^ 

The most important property to note is that the off- 
diagonal matrix element vanishes as gl"^™! as g (this 
behavior comes from the exchange term). Secondly, by 
exchanging m and n, it can be seen that 7i++ is a sym- 
metric matrix. Now for the second kind of off-diagonal 
element 



and pomiq) = {Pma(—<l))*- It is clear that pio ~ q, and 
P20 ~ a-iid the other matrix elements go with higher 
powers of q. Therefore one needs to keep only the com- 
ponents 10, 20, 01, and 02 of all the eigenvectors to 
find S{q) to order q"^. 

In order to compute these components of the eigen- 
vectors, we need to analyse the small-g structure of Ti in 
greater detail. In the following, we will set 6', = 7r/2 for 
simplicity. This choice makes TL real. 

First consider the diagonal matrix elements. 

7^(m0;m0;,) = e{m) - ^(0) + ^-^e-^^^^ 

- ! ^)V{s)e~^'/^L^{^^)MQS) (64) 

Let us recall that the differences of energies e(m) — e(0) 
are given in terms of the potential (eq(|2l|)) 

6(m) - e(0) = / ^„(.)e-^'V2(|^„^(,-^|2 _ |p„„(,-^|2) 

= (65) 

Consider 7i(10; 10, g). For small q the direct term 
contributes only in order q (for a long-range potential, 
v{q ^ 0) ~ 1/g), or order q^ (for a short-range potential, 
v{q — > 0) ~ const) ^ while the energy difference and the 
exchange term contribute to order q^ . Setting Jo ~ 1 
to get the order 1 contribution, we find that it vanishes 
for an arbitrary potential! The other diagonal matrix ele- 
ments mO for m 7^ 1 show no such behavior. This means 
that the naive magnetoexciton Oio has zero energy as 
q — + 0. This signals us that it must be related to the 
constraint. The detailed manner in which it is related 
will be seen shortly. 

As an aside we note that when one tries to construct 
CF-magnetoexciton wavefunctions in the standard man- 
ner by taking a IQH magnetoexciton, multiplyingJav the 
Laughlin- Jastrow factor, and projecting, one findsc3 that 
the IQH 10 magnetoexciton projects to zero as q — > 0. 
This seems to be directly related to what wc find here, 
but since it is difficult to go back to electronic wavefunc- 
tions from the Hamiltonian theory, it is not clear how to 
extend this relationship further. 

Let us now consider the off-diagonal elements of Ti.. 
There are two different types, those confined to the posi- 
tive energy subspace, and those that go between the pos- 
itive and negative energy subspaces. 



H(mO;On;g) = (-1) 



_u(g) 



^ /2 ( cQ\">-+n 



(67) 

Once again, this vanishes as (long-range poten- 

tial) or (short-range potential). 

Thus, as q ^ 0, the naive magnetoexcitons OmOj Oom 
become decoupled for m > 1, and therefore become the 
true eigenstates of Ti.. Oio and Oqi have a diagonal en- 
ergy of order q (long-range potential), and are coupled by 
matrix elements of order q, so they do not decouple. We 
will use these g — > naive magnetoexcitons as labels for 
the true eigenstates, denoting the positive energy ones by 
'i'm and the negative energy ones by ^'-m (for m > 1). In 
order to find the true eigenstates ^'q at small q, one can 
do perturbation theory with q as the small parameter. 

To see which of the true eigenstates we need to keep, 
consider the components of V&a. We have already seen 
that we need to restrict attention to the components 10, 
20, 01, and 02, due to the powers of q already present in 
p. Due to the small q structure of H it can be seen that 



|30 > +aq\2Q > +(3q\A0 > +-fq'^\m > 



(68) 



where |mO > represent the naive magnetoexciton states. 
Therefore the overlap < p{q)\'^§ >« g'^, and it is unnec- 
essary to compute it if one wants S{q) to order g*. Similar 
reasoning holds for Vl/m for any m > 2 and m < —2. 

To summarize, we see that to obtain S{q) to order 
it is sufficient to compute | < p{q)\'i'2 > P, since 'ifi will 
belong to the null (or unphysical) subspace. Also, the 
overlap needs to be computed only to an accuracy of q^. 

To make further progress, let us write down the matrix 
Ti. restricted to the lowest-energy four-dimensional sub- 
space spanned by |20 >, |10 >, |01 >, and |02 >, where 
all the diagonal matrix elements have been kept to order 
g^, and a long-range interaction has been used. 



Bk 
0(P) 



Bk 
Ak 
-Ak + Ck^ 
0(P) 



0(P) 
Ak - Cfc2 
-Ak 
Bk 



0{k^) 
0{P) 

Bk 
-E2o{k) 



(69) 



We have defined k = ql. The various coefRcients are 
(for the Coulomb interaction with the energy scale e'^ /el 
removed, v{q) — 2'!Tl/q) 



8 



A = 
B = 
C = 

£^20 — 



1 

3 
1 

4 



q+0{q^) 



(70) 



The central 2x2 submatrix consisting of the states 
10 and 01 can be seen to be very singular. Taken by 
itself, it possesses two zero eigenvalues (up to order q^) 
but only one eigenvector, and represents the constraint 
in its small-g form. The larger matrix also possesses two 
zero eigenvectors. This can be seen from the determinant 
of 7^4, which is 



det\HA\ = 2AE2Qk^{CE2o - B^) + k'^{B'^ 



^2 77i2 



The key identity 



=CS2o(0) 



(71) 



(72) 



makes this determinant zero to order fc"* inclusive, which 
means that there are two zero eigenvalues to order k^. 
Eq(^) seems accidental, and dependent on the potential 
used, but it can be shown from the form of the low-order 
Laguerre polynomials that this identity holds indepen- 
dent of the potential, and furthermore that 



B 



2^/2 



£^20(0) 



(73) 



for any potential. 

As mentioned before, the true eigenvectors can be 
found by straightforward perturbation theory. We first 
diagonalize the positive energy subspace by a unitary 
transformation (the Hamiltonian is hermitian in this sub- 
space). The eigenvectors can be used to construct the 
unitary transformation which needs to be correct only to 
order k 



1 



Bk 

'-' — I Bk_ 2 

E20 



(74) 



The resulting eigenvalues are E2 



E and Ak 



^20 I —e — 

B^k^ /E20 = Ak — Ck^, and the matrix in the new basis 
looks like 



U^HiU 



E 

Ak- Ck^ 

0{k^) -Ak + CP 
0(fc3) 0(fc2) 



0(fc2) 
Ak - Ck^ 
-Ak + CP 





(75) 



It is clear that to order q the nonzero positive eigenvec- 
tor is decoupled from the negative energy subspace. The 
central 2x2 matrix is degenerate, and has two zero eigen- 
values, but only one eigenvector. It represents the null 
subspace of Ti. in the new basis. In order to compute 
S{q) to order q"^ we need only the physical positive en- 
ergy eigenvector ^'2. It can be seen that ^'2 is decoupled 
from the null subspace. To the requisite order ^'2 is 



1 



o{k') o{e) 



(76) 



It can be easily verified that 5*2 is orthogonal to the null 
eigenvector, and this fact provides a shortcut to deriving 
the above results. Now taking the overlap with \p > and 
squaring, one easily finds 
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(77) 



This coefficient differs by a factor ojE.2 from that pre- 
dicted from Laughlipis wave functiorjj using the com- 
pressibility sum ruleB, namely S{q) = {ql)^/A. There- 
fore this calculation, while unambiguously showing the 
signature of an incompressible state, does not give the 
same result as the Laughlin state. This is a little sur- 
prising, since one expects the CF-HF state of one CF-LL 
filled to correspond to the Laughlin state. However, two 
points should be noted here. Firstly, it is very difficult to 
make the correspondence between the CF-wavefunction 
and the electronic one in the Hamiltonian formalism. It 
may be that the electronic wavefunction corresponding 
to one full CF-LL is incompressible, but not Laughlin's 
wavefunction. Secondly, this is an approximate calcula- 
tion, albeit one that is conserving, and a more sophis- 
ticated conserving approximation may well produce the 
same answer as Laughlin's wave function. 

Finally, one may extend this result to all the Laughlin 
fractions 2^+1 attaching 2s units of flux, with = 



2s 
2S-I-1 



to obtain 



S{q) = M' + ' 



(78) 



independent of s, which also differs from the result fai' 
the coefficient obtained from Laughlin's wavefunctionES, 
namely s/4. 

It should be emphasized that the results are indepen- 
dent of the form of the potential. The short-range po- 
tential produces a slightly different form of the matrix 
7^4 but one still obtains the same answer for the small 
q behavior of S{q). One expects this since the structure 
factor is a property of the ground state, and should not 
depend on the dynamics. 



B. The Structure Factor for the Principal Fractions 

ly =p/2p+ 1, p > 1 

Let us now turn to the more generic case oi p > 1. 
In this case we have p CF-LLs filled. We have already 
seen that only naive magnetoexcitons with A^f('ti2) 7^ 
Npirni) contribute to the density correlator, so we will 
only consider those below. Once again, we find that the 
diagonal matrix element Ti.{p,p — l;p,p — l;q), which is 



Ji^v{s)e-^'/^L,_,{^)L,{^)MQS) 



(79) 
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vanishes as q 
(eq@) 



0. Recall that the energy difference is where 



e{p)-eip-l) = f 



v{q)e 



p-1 
k=0 



B2 



~ D 

Vl2^2 — (^+- 



Ei)B2 
Ei)Bi 



D 



\ p—k 



{LI 



(80) 



The p, p— 1 diagonal matrix element vanishes as g ^ for 
an arbitrary potential by virtue of a Laguerre polynomial 
identity 



E 

fe=0 



k\ 



L.p^i{x)Lp{x) 



(93) 
(94) 

Let us call the basis vectors here ^f-, and ^Pq. Now 
it is easy to construct new basis vectors and '^q 

in the positive energy subspace that make it completely 
diagonal to order by doing first order perturbation 
theory in the off-diagonal elements. Once again it turns 
out that with this transformation, the positive energy 
subspace is decoupled from the negative energy and null 
subspaces to order q^ . The structure factor is then com- 
puted from the overlaps of the density with the physical 
(81) positive energy eigenvectors. 



This identity has been verified to high order p, but a 
general proof is not known to the author. 

Once again, this vanishing energy is related to the con- 
straint. Proceeding to the modes that are connected to 
this one by order q matrix elements, we find that there 
are two, Op+i^p-i and Op^p-2- This is the primary dif- 
ference between the case of generic p and p = 1, since 
in the latter case the second operator was missing. Now 
one proceeds as before, by first looking at the positive 
energy subspace 



/ EM Vi2{q) Bik 
n++,3 - Vi2{q) E2{q) B2k 
\ Bik b2k Eo{q) 



where 



Eo= n{p,p- l;p,p- l;q) 
Ei^nip+l,p-l;p+l,p-l;q) 
E2^ nip,p-2;p,p^2;q) 
Vi2^ nip+l,p-l;p,p-2;q) 
n(p,p - 1 ■,p + l,p - l;q) = Bik + ■ ■ ■ 
n{p,p-l ■p,p-2-q)^B2k + --- 



It should be noted that V12 does not vanish as g ^ 0. 
Next one carries out a unitary transformation that makes 
the upper left 2x2 submatrix diagonal. This unitary 
transformation is 



(82) 



(83) 
(84) 
(85) 
(86) 
(87) 




where 



E± = \{E^+E2± ^{Ei - E2Y + ) 
D = 



ViE+-E,] 



^12 



(89) 



(90) 
(91) 



After this transformation the Hamiltonian in the 3x3 
positive energy subspace assumes the form 




(92) 



S{q)^^^{\<p{q)\^+{q)> 



i?2^1~Vl2 -^2 \ ^ 

EiE2-V'i 



<p(g)|*_('?)>n 

1 I ■E'i.B2 — V12-B1 



EiE2-V:i 



(95) 



All the above algebraic steps are heavily dependent on 
the potential used, but nevertheless one expects that the 
final result should be potential independent. Once again, 
a formal proof is lacking, but this potential-independence 
has been verified for the Coulomb potential and the short 
range Gaussian potential for 2/5 and 3/7. By some mys- 
terious (to the author) properties of the Laguerre poly- 
nomials and their integrals, all the dependence on the 
functional form of the potential drops out. 

Let us illustrate this with the Gaussian potential 
v{q) = 27r(/*)^e^"'^ which is very convenient for ex- 
plicitly carrying out integrals, but which (for large a) 
also possesses a small parameter f3 = 1/(1 -I- a) which 
will be exploited later. For p — 2 the relevant quantities 
(without any approximation) are 



El = 


f - f /33 + _ 2^^5 


(96) 


E2 = 




(97) 


V12 = 


1 / 12 02 96 03 1 192 o4^ 
^\ 5 P 25P I25P ) 


(98) 


Bi = 




(99) 


B2 - 


/3^(l-f/3 + i/3') 


(100) 


expression appearing in the denominator of 


eq(|95|) 



EIE2 



*'l2 



P'{$P- 



2^^ n2 I 2^ 

-grA^ -I- -5! 



(101) 



It can easily be verified that both E1B2 — V12B1 and 
E2B1 — V12B2 are proportional to the above expression, 
giving ratios independent of (3. The structure factor de- 
pends on these ratios only, and is thus independent of 
/3. 

For p > 3 it becomes very difficult to compute -B;, 
and V12 exactly, even for the Gaussian potential. How- 
ever, if one assumes the potential independence of S{q)^ 
then one can use the following strategy: Since S{q) does 
not depend on it can be calculated in the limit of small 
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p. In this limit, one need only compute terms of lowest 
nontrivial order in /3 in Ei^ Bi, and Vi2. Keeping one 
more order as a consistency check one finds for generic p 



E2 = f3^c^p- 

Vi2= -Vp^ 

Bo = 



1 



1) 



l)-/3'''c'*(2p2-6p + 2) + - 
"T(/32c2-2/33c4(p-l) + . 
1/32(1 -/3c2(3p- 2) + •••) 



Vp^/3'(1-/3c2(3p-2)---) 



(102) 
(103) 
(104) 
(105) 
(106) 



Using the above expressions in eq(| 
result quoted in the introduction 



one finds the 



{qlYp^-ip"" 



3p 



p^ 



1 



(107) 



This expression, with its divergence as p ^ 
consistent with the resulted that for a problem equivalent 
to the I' = ^ problem, S{q) ~ g^log(g). As p ^ oo, 
the radius of convergence of the power series expansion 
of S{q) must go to zero (else the structu re fa ctor would 
diverge for a range of q according to eq(107)), and the 
p —^ oo limit does not commute with the q — > limit. 

Finally, it is interesting to compare the conserving pre- 
diction for the term of S{q) with the naive calculation 
using the preferred density in the HF state 



Snaive{q) =< F\Pp{q) Pp{-q)\^H F > 



(108) 



This latter number turns out to be p/8. Readn3 wrote 
the conserving density-density correlator as a sum of the 
naive correlator of the preferred density and a term rep- 
resenting the coupling of this preferred density to a gauge 
field. Making a similar decomposition here we see that 
the gauge-field contribution increases with p, and comes 
to dominate the small-g structure factor as p ^ oo. 



VI 



EXPLICIT RESULTS FOR i AND | 



To motivate the results presented in this section let, us 
consider the relation of the all-q Hamiltonian theoryO to 
the original electronic theory. The all-g theory is alge- 
braically consistent at all q and has the same algebra for 
the electronic density as the electronic theory. It is also 
consistent with the small-g theory that was developed 
earlier by canonically transforming the electronic Hamil- 
tonian. If one could find a sequence of canonical transfor- 
mations that exactly mapped the electronic Hamiltonian 
into this one, then one could be confident that the two 
problems are identical. In the absence of such a transfor- 
mation, all we can really be sure of is that the all-g theory 
has the same small-g behavior as the electronic problem. 
It is therefore of interest to compute the magnetoexci- 
ton dispersions for q not very small to see if it agrees 
with the numerical results for the electronic problem. In 
this section we will present results for the spin-polarized 
magnetoexcitons for i and | and show that the minima 
in the lowest energy magnetoexciton dispersion occur at 



the same values of ql as those found for the electronic 
problem. 

As shown in the previous section, for very small g, 
the naive magnetoexcitons decouple and become the true 
eigenstates oiH. As g increases they become increasingly 
coupled, and both level repulsion (between positive en- 
ergy states) and level attraction (between positive and 
negative energy states) manifest themselves, giving rise 
to complicated magnetoexciton dispersions. The matrix 
Ti. is infinite-dimensional, but in any numerical calcula- 
tion, only a finite matrix can be diagonalized. We saw 
that the lowest nontrivial result for S'(g) could be ob- 
tained by keeping at most a 6 x 6 matrix. As g increases 
more and more CF-LLs have to be kept to obtain accu- 
rate results. The accuracy of the results were checked 
by two different methods. Firstly, the number of CF-LLs 
kept was increased until the energy of the magnetoexci- 
ton was stable. Secondly, we know that at every g there 
should be two zero eigenvalues corresponding to the un- 
physical sector. For a finite number of CF-LLs kept, 
and g not too small, these would-be zero eigenvalues are 
not exactly zero, but much smaller than any other eigen- 
values. The number of CF-LLs kept in the calculation 
was increased until these null eigenvalues were at least 4 
orders of magnitude smaller than the smallest physical 
eigenvalue. 



0.16 p 

o.u 

0.12 

0.1 - 

0.08 - 

0.06 - 

0.04 - 

0.02 - 









FIG. 1. Spin-polarized magnetoexciton dispersions for 
u — ^ a.t X = 0.91 and X — I. All energies are in units of 



The two figures are for the long-range potential v{q) — 
22|_g-9 A /2 jT^j, various values of A. It is worth noting 
that though the depths of the minima are A-dependent, 
and different from what are found in numerical diago- 
nalizatioBS or from CF-wavefunctions, the positions are 
correctuTl. This indicates that the all-g theory does cap- 
ture some of the important physics of the electronic prob- 
lem even at fairly large g. Furthermore, the results are 
calculated in the lowest level conserving approximation. 
A more sophisticated conserving approximation incorpo- 
rating the screening of the potential may produce better 
results at large g. 
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FIG. 2. Spin-polarized magnetoexciton dispersions for 
1/ = I at A = and A = 0.5L All energies are in units of 



VII. EXTENTION OF THE FORMALISM TO 
NONUNIFORM DENSITY 

In this section we will see how the fractional quasipar- 
ticle charge emerges from the formalism, provided flux at- 
tachment is carried out for thfeinhomogeneous problem. 
The all-g Hamiltonian theorycS is based on the expecta- 
tion that the ground state is uniform. In the presence 
of a localized quasiparticle this is clearly not a good as- 
sumption. More generally, it would be useful to have a 
conserving formulation which was applicable to the case 
of a nonuniform ground state, since this would enable a 
treatment of edge^ or disorder. 

To construct such a formalism, let us go back to the 
defining equations for the guiding center coordinates of 
the electron and the pseudovortex eq(^,^. Let us as- 
sume that the ground state density is n(r) = n -I- in(r). 
According to the principles of flux attachment, the CFs 
should see an additional effective gauge field correspond- 
ing to this nonuniform density which satisfies 

V X (5a = 47r(5n(f) (109) 

This additional gauge field should be present in the CF- 
velocity operator 11, leading to the commutation relation 

[n,„,n,0] = ^(1 - c2 - 47r/25ri(rl)) (110) 

With this change, the quasiparticle charge can be com- 
puted by a very simple argument, iclejitical in spirit to 
the counting argument given by JairH. Consider a re- 
gion with large area A which contains one extra localized 
quasiparticle (in addition to p filled CF-LLs) in its inte- 
rior. The total charge in this region (after subtracting 
the background) is 

Q = -e<5n(g ^ 0) = -e(l + (5n7?((7 ^ 0)) (111) 

where the 1 counts the extra CF that has been intro- 
duced, and iufi is the response of the filled CF-LLs to 



the introduction of the extra CF. Before the extra CF 
was introduced, all the states in the CF-LLs 0, 1, • • • , p— 1 
were occupied. If one assumes that the final state is adi- 
abatically connected to the initial one, then the same 
states continue to be occupied. However, the total ef- 
fective flux through the area A has decreased in magni- 
tude by 47r(5n(g — > 0). This means that each CF-LL has 
25n{q 0) fewer states in the area A than it did before. 
In essence, the states in the regions of smaller effective 
field have a larger cyclotron radius, so they expand and 
push out the other states. This means the total charge 
in the area A is 

Q ^ -eSn{q ^ 0) = -e(l - 2pSn{q 0)) (112) 

which can be seen to give 

The inclusion of the varying gauge field in 11 is essen- 
tial for the success of this argument. Unfortunately, this 
modification of 11 has an undesirable effect: Using the 
modified 11 in the expressions for the guiding center co- 
ordinates Re (eq(||)) and Ry (eq(^) spoils their commu- 
tation relations (eqs([8[|9|,[l0|)). This is a serious problem, 
since this means that the density operator constructed 
from the new R^ does not obey the magnetic translation 
algebra, and thus cannot be identified with the physi- 
cal electronic density. The density also fails to commute 
with the constraint, making a conserving approximation 
infeasible. 

However, there is a way to modify Re and Ry so that 
the algebraic consistency of the theory is restored. Of the 
three commutation relations (eqs(p|,p|,p^) only the first 
two are essential to having a consistent theory. Satisfying 
eq(^ would mean that the density operator constructed 
from Re could be identified with the electronic density, 
while satisfying eq(^) would mean that this density com- 
mutes with the constraint. If eq(p^ is not satisfied it 
means that the constraint algebra does not close. New 
constraints would be generated by commuting old ones. 
This is not a problem for a conserving approximation, 
since by the Jacobi identity the electronic density opera- 
tor (and thus the Hamiltonian constructed from it) would 
commute with the new constraints as well. 

It turns out that operators Re, Ry which satisfy the 
conditions of the above paragraph to first order in ^ can 
be constructed. They are 

Rea — 1"a ~^ JT^-^a ~ c(l+c)^ I ^^{^ ^ t) X 



{Xe.Sn{r- ^X) + Sn{r^ (114) 

Rva ^ T(x^ c(l+c) ^ c^(l+c)'-* / ^^(c + i) X 



{XJn{^~ + - (115) 

where 
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Sn{r~ zX) 



£0/311/3 



zX) 



(116) 
(117) 



and Sn{q) is the Fourier transform of the ground state 
expectation value of the electronic density operator, and 
is a commuting number. The 11 operators appearing in 
the above formulas contain the spatially varying gauge 
field Sd. The preservation of the commutation relations 
eqs(^lj^) to first order in 6n/n can be verified in a straight- 
forward manner. This is therefore a good representation 
of the electronic problem for small deviations of density 
from a liquid state. Presumably eqs( 114 ,ll£) can be ex- 
tended to all orders in Sn/n while still maintaining the 
crucial commutation relations, eqs(|8[||). 

The modified density operator can in principle be used 
to calculate the detailed charge profile of a quasiparticle. 
One can create an extra CF in the p*'* CF-LL, and solve 
the CFHF equations with the above expression for the 
density operator. There is a self-consistent aspect to the 
formulation: The 6n(q) that enters the density operator 
has to ultimately be the expectation value of the same 
density operator in the ground state. 

The modifications to Re and i?„ to order 6n/n (other 
than the change in H) are not relevant for the charge of 
an isolated quasiparticle, but they are crucial if one wants 
to compute bosonic response functions for a nonuniform 
state in a conserving approximation. This issuej4fill be of 
particular importance in the edge state problemc3|-J"here 
one knows that there are physical gapless modesE3, and 
separating them from the unphysical modes can only be 
reliably carried out within a conserving approximation. 



VIII. DISCUSSION AND CONCLUSIONS 

The fractional quantum Hall statesB have long 
the peculiar status of being understood in principld3€, 
and yet being resistant to analytical calculation of their 
response functions. In the past few years much under- 
standing has been acheived on how to compute these re- 
sponse functions. This understanding has largely been 
based on the notion of field-theatetic; flux attachment via 
a Chern-Simons transformatiaM^I E3. Outstanding issues 
in the initial CS formulationsO were partially solved by 
Shankar and the author by adding variables and carrying 
out a sequence of canonical transformations rtoj-djecouple 
the high and low energy degrees of freedonJl3^E^. How- 
ever, the resulting theory, though it had many desirable 
features at small q, was unsatisfactory when viewed as 
a theory of ±he Lowest Landau level. Shankar posited a 
formulationES which was algebraically consistent at all q 
as a LLL theory, and had all the desirable features of the 
small-g formulation. i_. 

In the past, analysis of this ail-q formalismtJ was lim- 
ited to the HF approximation in which one used the pre- 
ferred density pp = pe — c^Pv The preferred density 
makes explicit many nonperturbative properties of the 
true density operator, such as the fractional quasiparti- 
cle charge and the order q^ matrix elements out of the 



ground state for incompressible states. The price paid is 
that the constraints, once they are used to construct pp, 
are ignored, and one is not sure how robust the results 
are under a proper implementation of the constraints. 

The most important objective acheived in this paper 
is that the nonperturbative features-that were introduced 
by hand into the preferred densityll3 are seen to emerge 
very naturally from the theory when constraints are con- 
sistently imposed and flux attachment is carried out for 
states of nonuniform density. 

The conserving approximation used is identical in 
spirit to the one used by ReadO to obtain response func- 
tions for the V = 1 boson problem based on the formalism 
of Pasquier and HaldanecJ. The case considered here is 
in some ways more complicated in that it involves Lan- 
dau level states instead of plane wave states, but also 
in other ways simpler, in that finite-dimensional matrix 
equations have to be sohied to obtain the response rather 
than integral equationscB. 

The order q^ matrix elements of the charge density 
operator emerge when unphysical states are discarded 
and only physical states are considered. The crux of the 
matter is that for = only the p—1 to p transition 
can produce an order q matrix element in violation of 
Kohn's theorem, and as g — > this magnetoexciton state 
becomes identical to the constraint to order q, and thus 
unphysical. 

Once the formalism is extended to include variations 
of densitj£.iB,the ground state, the fractional quasiparticle 
chargeHOS also emerges from the theory. 

Finally, a look at the spin-polarized magnetoexciton 
dispersions for q not too small shows that though one 
does not yet have a rigorous connection between the all-q 
theory and the original electronic theory, the all-q con- 
serving approximation manages to capture much of the 
important physics, even at large q. 

In summary, it appears that the all-g Hamiltonian for- 
mulation of the fractional quantum Hall problem coupled 
to a conserving approximation is well-suited to comput- 
ing bosonic response functions in the spin-polarized case. 
An analysis of the edge and disordered problems also 
seems possible in this approach using the extention to 
small density variations presented in Section VII , and 
will be pursued in future work. 
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X. APPENDIX I 

Mast of the material in this appendix appears in 
ref. |M| , but is presented here to make the discussion self- 
contained. In order to compute the sums over CF-LL 
indices of products of pe or p^ matrix elements, it turns 
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out to be convenient to express all the coordinates in 
terms of the guiding center and cyclotron coordinates of 
the CF, which are 



7?^ 



CF 



jri^z X n 



which have the commutation relations 



(118) 
(119) 



Using the completeness of the cyclotronic states |m > we 
can recast this sum as 



m2 > 



(135) 



We can then use the commutation relations of eq(120) to 
combine the exponentials and obtain 



X! X™ 1 ™ ( ^1 ) Pmni2 ( 9*2 ) 



CF,ay RcF,!^] 



^CF,a' -^CF,p 
RcF,ai RcF,f3 



■ I'' 





(120) 
(121) 
(122) 



We can now express the electron and pseudovortex guid- 
ing center coordinates in terms of the CF operators R'^p 



and R^p as 



Re 

Rv — R^ 



^CF 



CF 



cR 



CF 



CF 



(123) 
(124) 



< mi|e-*(T-+=92)-i?^F|^2 > (136) 

Finally we separate the exponentials in the reverse order 
to get the second useful identity 

Note that the cyclotron parts of pe and do not com- 
mute, even though the entire operators, including the 
CF-guiding center parts, do. 



The CF guiding center and cyclotron coordinates, in 
turn, can be expressed in terms of two sets of canoni- 
cal, commuting creation and annihilation operators 



r: 



CF,x 



iR. 



CF,y 



2 „ 



R. 



CF,x 



R^ 



R^ 

^CF,x 



iRcF,y — ^ 

^^CF,y ~ ^ 

iR^ — / 

^^CF,y ^ ' 



T^4 

1 — c 



^4 



where 



[ac,al] = [ag,al] 



(125) 
(126) 
(127) 
(128) 



(129) 
(130) 



We can abstractly set up the states of the CF-cyclotron 
creation and annihilation operators 



|n>=^|0> 



(131) 



In terms of these states the matrix elements of eq([l9|) 
and eq(pO|) are simply expressed as 



,Aq) = < mi\e 



< TOi e = 



|m2 > 
|m2 > 



(132) 
(133) 



In other words, only the cyclotron part of the electron 
and pseudovortex density enters the matrix elements de- 
fined above. An important result we will need is the 
following sum 



(gi)Pmm2(92) 



(134) 
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